Field correlations and the ultimate regime of turbulent convection 



(N 

O 

(N 



i : 

^ ■ 
Oh' 



> 
00 

cn 
rn 

o 



Mahendra K. Vcrma/'El Pankaj K. Mishra/ Ambrish Pandcy/ and Supriyo PauP 

^Department of Physics, Indian Institute of Technology, Kanpur 208016, India 
^Computational Fluid Dynamics Team, Centre for Development of Advanced Computing, Pune 4II 007, India 

(Dated: January 15, 2013) 

Using direct numerical simulations of Rayleigh-Benard convection (RBC) under free-slip boundary 
condition, we show that the normalized correlation function between the vertical velocity field 
and the temperature field, as well as the normalized viscous dissipation rate, scales as Ra~^'^^ 
for moderately large Rayleigh number Ra. This scaling accounts for the Nusselt number {Nu) 
exponent to be around 0.3 observed in experiments. Numerical simulations also reveal that the 
above normalized correlation functions are constants for the convection simulation under periodic 
boundary conditions. 
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Thermal convection, which is ubiquitous in engineer- 
ing and natural flows, exhibits various interesting phe- 
nomena, like instabilities, chaos, spatiotemporal pat- 
terns, and turbulence [J, Rayleigh-Benard convection 
(RBC) is an idealized model of thermal convection in 
which a fluid is placed between two horizontal conduct- 
ing plates, with the lower plate being hotter than the 
top one. The dynamics of the flow is governed by two 
non-dimensional parameters: the Rayleigh number Ra, 
a measure of the strength of the buoyancy force, and the 
Prandtl number Pr, a ratio of the kinematic viscosity 
and the thermal diffusivity. One of the most important 
theoretical and technological problems that remains un- 
solved in this field is nature of heat transport in RBC 
specially for very large Rayleigh number; this is the sub- 
ject matter of this paper. 

Experiments and numerical simulations reveal certain 
universal properties for the heat flux [ll-Q. Various ex- 
periments show that the Nusselt number, which is the 
ratio of the total (convective-|-conductive) heat flux to 
the conductive heat flux, scales as Ra^ with (3 around 
0.3 for moderately large Rayleigh numbers [7-14|. How- 
ever for very high Rayleigh number (called the "ultimate 
regime"), Kraichnan predicted that (3 = 1/2; several 
exp eriments found no evidence of the ultimate regime 0- 
llj, while some others claimed its existence [isl Il4l|. In 
this paper we show that the velocity-temperature cor- 
relation and the viscous dissipation rate vary with the 
Rayleigh number so as to yield the Nusselt number ex- 
ponent to be 0.3 for the intermediate Rayleigh numbers. 

One of the earliest attempt to understand Nusselt 
number scaling was by Kraichnan Q who derived us- 
ing the mixing-length theory that Nu cx Ra}/^ for large 
Pr, Nu oc [PrRaY/^ for small Pr, and Nu ~ 1 for 
very small Pr. For very large Rayleigh number or 
the ultimate regime, Kraichnan Q showed that Nu oc 
[i?a/(lni?a)]i/^ Malkus 3 argued that the Nu expo- 
nent of 1/3 is due to the property that the heat flux is 
independent of the cell height. By separating the dissi- 
pation rates in the bulk and the boundary layers, Gross- 



mann and Lohse Q showed that for the bulk dominated 
convective flows, Nu ^ {PrRaY^^ when A„ < Xe, but 
Nu ~ Ra}/^ when A„ > Xe- Here A„ and Xe are the 
widths of the viscous and thermal boundary layers, re- 
spectively. The parameter space of validity for the above 
scaling has been detailed in Grossmann and Lohse [Hj. 
Shraiman and Siggia [6| , Castaing et al. 0] , and Gioni et 
al. computed the Nu exponent and deduced it to be 
2/7 due to the boundar y la yers. Using scaling arguments 
Verzicco and Camussi [l5j claimed that Nu ^ {Ra)^/'^ 
for small Pr {Pr < 1). 

Many experiments have been performed to test the 
above scaling, yet they have not been able to resolve the 
scaling exponents completely. Laboratory experiments 
on typical fluids, water, helium gas, and mercury, yield 
Nusselt number exponents from 0.26 to 0.31 for Ra up 
to 10^^ 0-[i3- Gioni et al. d and Glazier et al. ^ 
used mercury, while Gastaing et al. 0] and Niemela et 
al. used helium gas as their experimental fluid. Gioni 
et al. @ also performed experiments on water for com- 
parison. Ahlers et al. [l^ employed pressurized mixture 
of gas in their experiment. Fig. [T] illustrates the plots of 
reduced Nusselt number Nu/ (PrRa)'^''^'^ vs. PrRa com- 
puted in earlier RBG experiments and numerical Simula 
Ghavanne et al. 



14| reported 



tions. Ghavanne et al. |13j and Roche et al. 
existence of the ultimate regime in their experiment on 
helium. Roche et al. used nonsmooth surfaces to cancel 
the thickness variation of the viscous sublayer. Ahlers et 
al. [l2| found simultaneous existence of multiple Nusselt 
number exponents (0.17, 0.25, and 0.36) for large Ra's, 
possibly due to multiple coexisting attractors. 

The Nusselt number for RBG is also investigated using 
direct numerical simulation. Verzicco and Gamussi [l5| 
and Kerr and Herring [3| showed that Nu Ra^^^ for 
small Pr, while Silano et al. 17 1, Stevens et al. [l8| . and 



Verzicco and Sreenivasan [l9| reported the exponent be 
around 1/3 for very high Rayleigh number. Kerr (20j and 
Kerr and Herring [3| found the exponent to be around 
0.28 2/7 for larger Pr. Fig. [1] illustrates some of these 
results, as well as our numerical results for Pr = 0.2 and 
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FIG. 1. Reduced Nusselt number {Nu/ {PrRaf'^'') versus 
PrRa, the product of Prandtl number (Pr) and Rayleigh 
number {Ra). Experimental data: Cioni et al. for wa- 
ter (black v)- Numerical data: Verzicco and Camussi [l^ 
for Pr = 0.7 (green x); Kerr [2^ for Pr — 0.7 (magenta o); 
Stevens et al. for Pr — 0.7 (black •); Verzicco and Sreeni- 
vasan [l^ for Pr = 0.7 (grey hexagons); our simulation data 
for Pr = 6.8 (red □) and for Pr = 0.2 (blue *). 



tance d and with a temperature difference of A, are 

Stu + (u • V)u = -— + age*! + :/V^u, (1) 
Po 

dtB+{u-V)e ^^u^ + kV^O, (2) 

where 9 and a are, respectively, the temperature and 
pressure fluctuations from the steady conduction state 
(T = Tc + 0, with Tc as the conduction temperature pro- 
file), z is the buoyancy direction, po is the mean density of 
fluid, g is the acceleration due to gravity, and a, v, and k 
arc the thermal heat expansion coefficient, the kinematic 
viscosity, and the thermal diffusivity of fiuid, respectively. 
The two most important nondimensional parameters of 
RBC arc the Rayleigh number Ra — ag/S.(P /vk and the 
Prandtl number Pr ~ v j k. The Nusselt number can be 
expressed as 



Kl^/d+{u^T) 



= 1 



(3) 



where the nondimensionalized vertical velocity and tem- 
perature fields are u'^ = Uzd/n and 9' = 9/ A. 

Eq. ([21) has three competing terms: (u- V)6', ^Uz, and 
k\/^6. When the diffusion term kS/^9 of Eq. ^ is much 
smaller than the other two, we obtain 



Ul0l 



{u,)lA 



d d 
In the momentum equation [Eq. ^] 



A. 



d 



Ul ~ \J agAd, 



(4) 



(5) 



6.8 (to be described later). Kerr [2l| also studied the 
energy budget in RBC by computing the mean square 
velocity and dissipation rates as a function of Rayleigh 
number for various aspect ratios and Prandtl numbers. 

For bulk-dominated convective flows, a simple argu- 
ment would predict that Nu oc (u^ 

perature fluctuation). This scaling argument assumes 
that ^ Ra^f^ and {9^^^ - const, and it ig- 

nores the correlation between the vertical velocity field 
and the temperature field. Similarly, one of the exact re- 
lations Nu — 1 ^ /{PrRa) Q yields Nu ~ Ra^^'^ if we implying that the Peclet number Pe = ULd/n scales as 
replace e" ^ U^/L ignoring correlations {Ul is the large Pe « PrRa. These scaling relations are applicable 
scale velocity, which is the free-fall velocity). A major 
difficulty in this field is how to reconcile the Nu ^ Ra^'^ 
scaling observed in the experiments to the Nu ^ Ra^^^ 
predicted for the ultimate regime. In this paper, we will 
explicitly compute the velocity-temperature correlation 
function and show that it scales with Ra in such a way 
that Nu ~ Ra'^'^ at moderate Rayleigh numbers. 

Properties of convective flow depends quite critically 
on the boundary layers (2^ . However, the Nu scaling 
appears to be somewhat insensitive to the presence of 
boun dary layers [2^ and change of the boundary condi- 
tions [l9|. Motivated by the above observations, we at- 
tempt to compute the scaling of the Nusselt and Peclet 
numbers in turbulent regime in terms of bulk quantities, 
specially the large scale velocity Ul and the large scale 
temperature 9l {— \/ {9'^)). The combined effects of the 
bulk and boundary layers would be studied later with 
more refined simulations and theoretical arguments. 

The RBC equations under Boussinesq approximation 
for a fluid placed between two plates, separated by dis- 



when {uz)LA/d ^ kV^9 or PrRa ^ 1. These results are 
consistent with the predictions of Kraichnan Q , Gross- 
mann and Lohse [H , and Silano et al. [l^ . 

Note however that for PrRa <C 1, the diffusion term 
dominates the nonlinear term, and it is balanced by 
UzA/d. Therefore, 9l « C/l^A/k and Pe « PrRa. 
Under this condition, Ul oc 9l, hence Nu = 1-1- 
c{PrRa)^ with c as a constant. Convective turbulence 
with PrRa <ti 1 limit is very rarely observed in terres- 
trial experiments or astrophysical observations. There- 
fore, in the present paper we limit ourselves to PrRa 3> 1 
regime. In Fig. [5] we plot the reduced Peclet number 
Pe/ {PrRa)'^'^ computed in various RBC experiments 
and numerical simulations; these results are in general 
agreement with Pe ^ {PrRaY^^ scaling The data 
of Figs.[l]and[2]appear to be compactified reasonably well 
with PrRa, hence we use PrRa as a scaling variable. 

As discussed earfier, a naive replacement of (u^,^') of 

Eq. (01) with (M;2^)l/2^5,/2^1/2 yj^yg jy^ _ ^/pFR^, which 

is not observed in the experiments and simulations for 
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plot them in Fig. [3l We observe that 

C^^PrRa) « (5.6 ± 1. 30) (Pri?a) 



-0.22±0.017 



(7) 



for Ra ~ 10®. The above correlation is due to the inter- 
actions between the bulk flow and the boundary layers. 
However its power law behaviour is interesting which im- 
mediately yields Nu ^ Ra'^-^'^, consistent with the Nu 
scaling observed in experiments and numerical simula- 
tions. Also, a combination of Eqs. (jS]) and (O provides 



Nu-1 = 0M6{C''\PrRa)}{PrRa) 



1/2 



(8) 



FIG. 2. Reduced Peclet number {Pe/{PrRa)° '^) versus 
PrRa. Black v represent the experimental data for helium 
(Pr ~ 0.7) by Niemela et al. Red □ and blue * represent 
our simulation data for Pr — 6.8 and Pr = 0.2, respectively. 
The inset shows the constancy of the large scale temperature 
field {Ol/A) with PrRa for Pr = 6.8 (red □) and Pr = 0.2 
(blue *). 



moderately large Rayleigh numbers (see Fig. [T]). To ac- 
count for the velocity-temperature correlation, we rewrite 
the Nusselt number as 

iV^ - 1 = «0') = C-'{PrRa){uf)\/'{9'')\f\ (6) 



where C'^iPrRa) = 



ized correlation function between the vertical velocity 
and temperature, with (}y and ()t representing the vol- 
ume and the temporal averages, respectively. Note that 
C^^ {PrRa) < 1 as a consequence of Cauchy-Schwarz 
inequality. We perform numerical simulations to com- 
pute the above correlation, along with the Peclet number, 
Nusselt number, and 6l, for various Pr and Ra values. 

We numerically solve Eqs. using pseudo-spectral 

method on grid-size ranging from 128^ to 512"^. We ap- 
ply free-slip boundary conditions for the velocity field 
and isothermal boundary condition for the temperature 
field. We use Runge-Kutta fourth-order scheme for time 
advancement, and compute the relevant quantities in the 
steady state. For further details of the simulation, refer 
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In Figs. [T] and [5] we plot the computed reduced 
Nusselt and Peclet numbers vs. PrRa for Pr = 6.8 
(red squares) and 0.2 (blue circles). The best fits for 
our data are Nu ^ (0.27 ± 0.04)(Pri?a)°-27±o oi ^nd 
Pe = (0.26 ± 0.04) (Pri?a)°-49±° °i for Pr = 6.8, and 
Nu = (0.39 ± 0.02)(PrPa)"-27±° °i and Pe = (1.04 ± 
0.20)(Pri?a)°-^3±o.O2 f^^, ^ q q2 ^he inset of Fig. [2] 

demonstrates constancy of 0l/A with PrRa. These re- 
sults are consistent with the earlier experimental and nu- 
merical observations, as shown in the figures. 

We compute the normalized correlation function 
C^^ {PrRa) for Pr ~ 6.8 and various Ra values, and 
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FIG. 3. Normalized correlation function C^^ {PrRa) 

{{<e')v/{{u'^)\l'^{0''^)v'^))t versus PrRa. Function (5.6 ± 
1.30) (Pri?a) 
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We can also deduce another important field correlation 
using one of the exact relations for the RBC flow Q: 



Nu~\^ 



v^Ra 



{Pef_ 
PrRa 



C'^{PrRa), (9) 



where the dissipation rate of the kinetic energy e" = 
{Ul/d)C^" {PrRa), with the function {PrRa), named 
as "normalized viscous dissipate rate", playing similar 
role as C^PrRa). Since Pe = a{PrRa)^/'^, C^"{PrRa) 
must scale as s{PrRa)~^ '^'^ similar to C^^ {PrRa) (Here 
a, s are constants). Thus 



Nu-1^ a^{PrRa)^/^s{PrRa)-"-^^ = a^ s{PrRa)°'^'^ 

(10) 

that immediately yields s « 15. 

The normalized correlation function C"^^ {PrRa) and 
the normalized viscous dissipation rate C"^ {PrRa) de- 
crease with increasing Ra since the fields tend to be- 
come more and more turbulent. A question is whether 
these correlations would continue to decrease and they 
will vanish as Ra — > oo, or they would saturate to some 
asymptotic values. For very large Rayleigh numbers, we 
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expect the turbulence to be fully developed, and, ac- 
cording to the Kolmogorov's theory for fully-developed 
turbulence, e" - (C/£/i) or {PrRa) ~ const. This 
would correspond to the "ultimate regime" proposed by 
Kraichnan [3[. Therefore, we conjecture that the func- 
tions {PrRa) or C"^ {PrRa) would be as shown in 
Fig. |4l a power law for intermediate Ra, and constant 
for very large Ra (ultimate regime). 



I 
1. 




FIG. 4. The purple circles with extended chained line depict 
constancy of the normalized correlation function C^^ {PrRa) 
for the periodic box {Pr — 1). Assuming that the convec- 
tive turbulence becomes fully-developed for very large Ra, we 
conjecture that the normalized correlation function would be- 
come a constant in the ultimate regime. The second curve is 
reproduction of the correlation function of Fig. 3 along with 
extended Ra"*^'^^ for moderately large Ra, and then a con- 
stant for the ultimate regime after some transitional Rayleigh 
number (here 10^^). 

The constant values of the C^iPrRa) and {PrRa) 
in the predicted ultimate regime can be computed given 
the transitional Rayleigh number beyond which the ul- 
timate regime is expected. From the present set of 
experiments 0413' '^^^ conclude that the ultimate 
regime possibly starts beyond Ra ~ 10^'^[l3|- With this 
Rayleigh number, the asymptotic values of C^^ {PrRa) 
and ^^^{PrRa) would be around 0.00066 and 0.0018, 
respectively. Using these values, we conjecture that the 
Nussclt number in the ultimate regime would scale as 
Nu w 3.1 X lQ-'{PrRa)^/'^. 

It has been conjectured that turbulent convection of 
the ultimate regime is approximately represented by con- 
vection with periodic boundary conditions at significantly 
lower Rayleigh numbers 27[. To strengthen our theo- 
retical arguments presented earlier, we carried out RBC 
simulations for the periodic box geometry of size (27r)^ 
using a pseudospectral code [25, 26|. We computed the 
normalized correlation function C"^^ {PrRa) for Pr = 1 
and high Rayleigh numbers Ra — 10'*-10*, for which 
the flow is fully turbulent. We observe that for all the 
runs C^^ {PrRa) ~ 0.75 as shown in Fig. HI as well 



as C'"{PrRa) « 0.48. Also, Pe « 5.7Ra°-^°^'° °'^, 
Nu w 23.7i?a°-46±° °4, and 9/ A w 4.21 ± 0.21. More- 
over, the flow is highly anisotropic with (2|uzp/(|uj,p + 
« 3.26 ± 0.43. Thus, our prediction for the uhi- 
mate regime is consistent with the convective flow with 
periodic boundary conditions. 

In summary, we relate the Nusselt number scaling 
to the normalized velocity-temperature correlation func- 
tion, as well as the normalized viscous dissipation rate. 
We show that these functions scale with the Rayleigh 
number as i?a~°'^^ for intermediate i?a's, which leads 
to Nu ^ Ra'^'^ observed in experiments. For very large 
Rayleigh numbers, the flow is expected to be fully tur- 
bulent. This observation leads us to conjecture that 
the normalized velocity-temperature correlation function 
could become a constant for very large Ra, thus yielding 
Nu '--^ Ra^^"^, as predicted by Kraichnan Q. Convec- 
tion simulations for a periodic box geometry also exhibit 
constant values for the normalized correlation functions, 
as well as Nu ~ Ra^^^, consistent with our conjecture. 
Thus, the field correlations provide valuable insights into 
the scaling of the Nusselt number. 

Numerical simulation of ultimate regime is beyond the 
capabilities of present computers, but gap is expected to 
be bridged soon, which would settle this long-standing 
problem. Also, attempts to compute the asymptotic val- 
ues of the correlations functions in the similar lines as 
the theoretical computations of Whitehead and Doering, 
and Constantin and Doering [28| would be very valuable. 
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MKV, and a BRNS grant. 



[1] 
[2] 



[3] 
[4] 



lmkv@iitk. acini 

E. D. Siggia, Ann. Rev. Fluid Mech. 26, 137 (1994). 
G. Ahlers, S. Grossmann, and D. Lohse, Rev. Mod. Phys. 
81, 503 (2009); D. Lohse and K. Q. Xia, Ann. Rev. Fluid 
Mech. 42, 335 (2010). 

R. Kraichnan, Phys. Fluids 5, 1374 (1962). 
M. V. R. Malkus, Malkus, Proc. R. Soc. Lond. A 225, 
185 (1954); M. V. R. Malkus, Proc. R. Soc. Lond. A 225, 
196 (1954). 

S. Grossmann and D. Lohse, J. Fluid Mech. 407, 27 
(2000). 

B. I. Shraiman and E. D. Siggia, Phys. Rev. A 42, 3650 
(1990). 

B. Castaing et al, J. Fluid Mech. 204, 1 (1989). 

S. Cioni, S. Ciliberto, and J. Sommeria, J. Fluid Mech. 

335, 111 (1997). 

J. J. Niemela et al. Nature 404, 837 (2000). 



5 



[10] J. Glazier, T. Segawa, A. Naert, and M. Sano, Nature 

398, 307 (1999). 
[11] P. Urban, V. Musilova, and L. Skrbek, Phys. Rev. Lett. 

107 (2011). 

[12] D. Funfschilling, E. Bodenschatz, and G. Ahlers, Phys. 

Rev. Lett. 103, 014503 (2009). 
[13] X. Chavanne et al, Phys. Fluids 13, 1300 (2001). 
[14] P. E. Roche et al. Phys. Rev. E 63, 045303(R) (2001). 
[15] R. Verzicco and R. Camussi, J. Fluid Mech. 383, 55 

(1999). 

[16] R. M. Kerr and J. R. Herring, J. Fluid Mech. 419, 325 
(2008). 

[17] G. Silano, K. R. Sreenivasan, and R. Verzicco, J. Fluid 

Mech. 662, 409 (2010). 
[18] R. Stevens, R. Verzicco, and D. Lohse, J. Fluid Mech. 

643, 495 (2010); R. Stevens, and D. Lohse, R. Verzicco, 

arXiv;1102.2307 (2011). 



[19] R. Verzicco and K. Sreenivasan, J. Fluid Mech. 595, 203 
(2008). 

[20] R. Kerr, J. Fluid Mech. 310, 139 (1996). 

[21] R. M. Kerr, Phys. Rev. Lett. 87, 244502 (2001). 

[22] M. van Reeuwijk, H. J. J. Jonker, and K. Hanjalic, Phys. 

Rev. E 77, 036311 (2008); Phys. Rev. E 77, 036312 

(2008); 

[23] R. Verzicco, EPL 35, 133 (2003). 
[24] J. J. Niemela, et al. J. Fluid Mech. 449, 169 (2001). 
[25] P. K. Mishra and M. K. Verma, Phys. Rev. E 81, 056316 
(2010). 

[26] M. K. Verma, Arxiv:1102.251 [physics.flu-dyn] (2011). 
[27] D. Lohse and F. Toschi, Phys. Rev. Lett. 90, 034502 

(2003); E. Calzavarini et al, Phys. Fluids 17, 055107 

(2005). 

[28] J. P. Whitehead and C. R. Doering, Phys. Rev. Lett. 
106, 244501 (2011); P. Constantin and C. R. Doering, J. 
Stat. Phys. 94, 159 (1999). 



